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Dielectric Barrier Discharge (DBD) plasma actuators with serpentine shaped 
electrodes cause a change in the operational behavior of the flow control mechanism 
relative to the standard ‘linear’ actuator due to the introduction of fully three 
dimensional vortices. A parametric study is performed on a partially separated 
transitional flow over an SD7003 airfoil. This study examines the change in behavior 
as a function of geometric parameters, in particular the geometric amplitude of the 
serpentine geometry as it is flattened out, reverting back to the standard linear 
actuator. Continuous operation of the actuator does not fully remove the 
separation, but it does modify the properties of the separated boundary layer. 
Separation is eliminated when the actuator is operated with a duty cycle. 


Nomenclature 
c = Chord length 
Cp = Coefficient of drag 
Cr = Coefficient of skin friction 
CL = Coefficient of lift 
CM = Coefficient of moment 
CG. = Coefficient of pressure 
D. = Nondimensional plasma force magnitude 
Ẹ = Energy 
Eo = Characteristic electrical field 
e. = Electron charge 
éG N = Body fitted coordinates 


i = Body force in the 1 direction 
Nondimensional forcing 
Inviscid flux vectors 

Viscid flux vectors 

Ratio of specific heats 
Jacobian 

= Mach number 

= Pressure 

= Density 


ee 
Qe 
= 
| | 


zas 


a Q 


OODD 


U, V, W 
u, V, W 


= Characteristic electron density 


Heat Flux 

Vector of dependent flow variables 
Contravariant velocity components 
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Re = Reynolds number 


S = Source term vector 

St = Strouhal number 

T = Nondimensional time 

uf = Shear stress 

X,Y, Z = Cartesian coordinates 
Subscripts 

5 = Freestream reference velocity 


= Partial derivatives with respect to the x;" direction 


xi 


I. Introduction 


Separated flows and the methods used to reattach them are two highly documented issues in flow 
control research’. The separated flows that can exist over airfoils have a number of performance reducing 
properties, including increased drag, stronger vortex shedding and reduced ratio of lift to drag, among 
other effects. Reattachment of these flows at low Reynolds number is of vital importance to many Micro 
Air Vehicle (MAV) applications. For these types of applications, the flow often has some level of laminar 
behavior, which is important for control of the flow, as transitioning these flows can lead to reattachment. 

DBD plasma actuators have a number of beneficial properties that are unseen in other flow control 
devices’”. The devices themselves have no moving parts, and respond on the time scale of an electrical 
device, far faster than the fluid itself. Plasma actuators can also be flush mounted such that they have no 
impact on the flow field when not in use. Furthermore, DBD plasma actuators have shown themselves to 
be effective for boundary layer control*® and for flow reattachment at low Reynolds number®'', 
especially when used in a pulsed manner, where unstable frequencies can be excited, causing the flow to 
transition and then reattach. 

Plasma actuators with novel geometries, such as the serpentine geometry have shown themselves to 
be very useful in a number of numerical studies''’* as well as limited experimental work’”'’. These 
actuators are able to introduce streamwise, as well as spanwise and surface normal vorticity to the flow, 
and promote increased mixing of the flow. 

This study aims to generalize the geometry of the plasma actuator electrodes and explore the flow 
physics as this geometry is modified. In this generalization, the standard actuator is considered to be a 
special case of the serpentine geometry (with infinite wavelength and zero amplitude) where no 
streamwise or surface normal vorticity is added to the flow by the actuator except for at the spanwise 
edges of the electrodes. While published results show an interesting difference in the physics when the 
electrode geometries are distinctly different’’, the numerical simulations performed for this work hope to 
fill the knowledge gap as the actuator is smoothly transitioned from the serpentine geometry to the linear 
geometry. 


II. Serpentine Geometry Plasma Actuation 


DBD plasma actuators are small, electrically powered devices composed of two asymmetrically placed 
electrodes, one exposed and one encapsulated in a dielectric material. A high electrical potential is 
applied between the two electrodes. The electrical field in this region causes gas on the surface near the 
two electrodes to become weakly ionized’. The ions and electrons in the flow are imparted momentum 
via the electric field, and through collisional processes, momentum is imparted to the remaining neutral 
fluid particles. 

Novel plasma actuator geometries, such as the finger geometry of Satchzman et al." or the plasma 
synthetic jets developed by Santhanakrishnan and Jacob'® have been used for various flow control 
applications. The serpentine plasma actuator developed by Roy and Wang'~”” has its own benefits in that 


2 


it is able to induce spanwise and streamwise vorticity 
(through the yx f term in the vorticity equation), 


while still maintaining relatively simple electrode 
geometry that can be extended indefinitely to any 
distance (Figure 1). This geometry is also significant 
in that it is able to define a characteristic spanwise 
wave length to the flow control, which is important 
with respect to exciting certain unstable modes in 
order to transition the flow. 


III. Governing Equations and Numerical Models 


This study was performed using the Implicit 
Large Eddy Simulation (ILES) method’*”’ using the 
implicit Beam and Warming algorithm’ for time 
integration. The unsteady, three dimensional, 
compressible, unfiltered form of the Navier-Stokes 
equations is solved, and then a Pade type filter such as those developed by Lele” is applied. The 
equations are solved in a generalized body-fitted coordinate system. The equations are solved in the 
conservation form, where 


<(59)+ <(r- —F,) +—(6 Gy) + —(H —H) =D S 
Ot QJ + ag Re ”) ôn Re ”) 0 Re) ce 


Q is the vector of dependent variables. The terms F, G, and H are the inviscid flux vectors, and F,, G,, 
and H, are the viscid flux vectors. Furthermore t represents time, and €, n, and C represents the body fitted 
coordinate system. The plasma forcing is included via the source term S, explained in more depth below. 
While this is an incompressible flow, because a compressible flow solver is used, a Mach number of 0.1 
is used in order to minimize compressibility effects while maintaining numerical stability. There terms are 
defined as the dependent variables as 
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Figure 1. Serpentine geometry from the side 
(a) and from above (b). From Durscher and 
Roy, 2011, used with permission. 


the inviscid flux vectors are 
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the viscid flux vectors are 


0 0 0 
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p represent the pressure. p represents density. T represents temperature, and finally, E represents 
energy. All of these variables have been nondimensionalized by their characteristic values, except for 
pressure, which has been nondimensionalized by the dynamic head of the characteristic variables. The 
length scales have all been nondimensionalized by the chord length of the airfoil. 

The plasma body force comes from the 
Townsend ionization of the air between two 
asymmetrically placed electrodes, separated by a 
dielectric material (Figure 2). The magnitude of this 
forcing is controlled by the nondimensional 
parameter D., being the ratio of twice the maximum 
possible plasma work to the dynamic head as shown 
below. A reduced order model for the plasma was 


— used. While this model is simplest of the available 
L models, the primary purpose of the plasma in these 
Figure 2. Schematic of a DBD plasma actuator simulations is not an in depth study of the plasma 
and the model of the body force. kinetics, but how disturbances due to plasma 

actuation affects the flow overall. Additionally, recent 
results!? have shown that there is relative independence of the plasma model to the end result under 


certain conditions, similar to what is used in this simulation in this simulation. The plasma source is 
formally described for these simulations as 
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The computations for this study were done on a cylindrical O-mesh of 390x242x97 (angular x 
radial x spanwise) points representing an SD7003 airfoil with a unit chord length and extended 0.2 chord 
lengths in the spanwise direction. The mesh extends radially outwards 100 chord lengths in order to 
provide a good free stream condition. The mesh is also stretched as it goes to the free stream in order to 
prevent reflections off of the free stream boundary. Periodic boundary conditions were applied in the 
spanwise and angular directions. The baseline flow case was also on a comparable mesh with 
650x390x151 points representing the same domain in order to check for grid independence of the 
solution. 





Figure 3. Geometries of the serpentine force distributions. From left to right - Full, three quarters, half, 
and quarter amplitudes of the serpentine geometry. 


Force Geometry Amplitude (r) Wavelength (p) 
Full 0.02500 0.1 
Three quarters 0.01875 0.1 
Half 0.01250 0.1 
Quarter 0.00625 0.1 


Table 1. Details of the plasma actuator geometries. 


Case Force Geometry Pulsing Freq. (St) D. 

Baseline None n/a n/a 
Baseline — fine None n/a n/a 
Full — Continuous Full 5.0 3:21 
Three quarters — Continuous Three quarters 5.0 3.27 
Half — Continuous Half 5.0 3.27 
Quarter — Continuous Quarter 5.0 3.27 
Full — Pulsed Full n/a 3-21 
Three quarters — Pulsed Three quarters n/a 3.27 
Half — Pulsed Half n/a 3.21 
Quarter — Pulsed Quarter n/a 3.27 


Table 2. Description of the plasma body force for the cases run 


A baseline flow for the SD7003 airfoil was established after running for 100,000 (10 non-dimensional 
units of time, t, where t=u,,/c) times steps (At=0.0001) to remove any transient irregularities. This flow 
was used as an initial condition for all of the cases tested. All of the cases were run for at least 5t in order 
to reach an equilibrium state where the transient effects could be considered negligible. Statistics were 
then collected over a period of 3.2 t (equivalent to 16 cycles of plasma forcing). 


IV. Results and Discussion 

A. Baseline Case 

A baseline case was run in order to establish a control case from which the others can be compared to. 
The solution to flow over an SD7003 airfoil at a Reynolds number of 40,000 with a 4° angle of attack 
shows a separated shear layer over the upper surface of the airfoil. This shear layer is laminar in nature, 
but roughly two thirds down the chord length the flow transitions and reattaches to the surface of the 
airfoil. The flow over the last two thirds of the airfoil continues to transition, and once the flow reaches 
the trailing edge of the airfoil, it is fully turbulent (Figure 4a). 

A region of low momentum exists in the separation bubble. This bubble is largely steady in nature, 
except for near the transition point (Figure 4b). This bubble is largely steady in nature, except for near the 
transition point (Figure 4b). 


This flow also shows a quasi-periodic shedding of 
vortical structures starting within the laminar shear 
layer with a frequency around St=5, indicating that 
some sort of instability exists near that frequency. 
Previous'’”'' have taken advantage of this instability to 
remove the flow separation and those studies were able 
to produce an attached flow. 

In addition to the structure of the flow, the size and 
shape of the separation bubble match sufficiently well 
with previously reported experimental and numerical 
data ^. 


Vel. Mag 


B. Continual actuation with serpentine geometry — 
flow and surface features 

Running the serpentine geometry with continuous 
actuation allows for continual streamwise vorticity to be 
added to the flow as a function of the geometry, 
isolating the effects of streamwise vorticity insertion 
from those of any pulsed flow stability, which are e 
addressed below in Sec. 4C. A 

As the amplitude of the serpentine geometry and the 
amount of streamwise vorticity is increased, the 
separated shear layer begins to show waviness to it, as 








the streamwise vorticity lifts it in certain areas, and 
depresses it in others (Figure 5). This waviness shows 
that the actuator is having an effect on the boundary 


Figure 4. Instantaneous velocity magnitude 
(a), time and spanwise averaged velocity 
magnitude (b), and streamlines of the time 


layer, even though the separated region is not attaching. 9d spanwise averaged flow (c). 


Previous studies” have shown this type of continuous serpentine actuation to be effective for 
separation control using more concentrated plasma forcing (over a smaller volume, but with the same 
total amount of force). However, in the present study, the shear layer does not break down, and separation 
control is not achieved. Because the forcing here is less compact than the forcing used in that study, a 
significant portion of it is inserted outside of the boundary layer. As much of the force is inserted outside 
the boundary layer, less of the total force is available to perturb the boundary layer itself and is used to 
perturb the free stream, which has a lesser effect on the flow, though enough of the body force is still 
present in order to perturb the flow. 

The time averaged flow field confirms that the flow is still separated, and further shows that the 
transitional effects are not strongly impacted (Figure 6). Additionally, the primary effect seems to be 
slight forward movement of the transition of the shear layer as the geometry of the serpentine actuator 
becomes more pronounced. Even still, the thickness of the wake and the velocity deficit seem to be nearly 
constant. The lack of streamwise vortical effects via the Q-criterion” indicates that the actuator is not 
causing streamwise rotational effects. Instead, it is introducing more shear stress than rotational stress in 
the streamwise direction. 

The time averaged surface pressure shows that a plateau of surface pressure exists for the 
uncontrolled case. This is typical of recirculating, separated flows. As the continuous plasma flow control 
is applied, the shape of this plateau changes, shortening in length, but slightly increasing in height as the 
serpentine geometry becomes fuller (Figure 7a). This change in the surface pressure profile suggests that 
the laminar boundary layer is thickening (causing a larger pressure drop), but transitioning slightly 
quicker. 





Figure 5. Velocity magnitude (at z=0.1 the pinching point, left) and Q-criterion (Q=100, right) for 
the cases of continual serpentine plasma actuation, varying the geometry of the actuator (from top 
to bottom - full, three quarters, half and one quarter). 


Subtracting the spanwise average from the time averaged of the surface pressure coefficient, one can 
see that there is a certain noisy periodicity to the surface pressure with the roughly the same spanwise 
wave number as that of the plasma actuator even after 100 samples (Figure 7b-e). This is another 
indicator that while the separation is not removed using continual actuation, there are subtle effects of the 
actuation on the flow. Using a stronger actuator would amplify these effects, and eventually control could 
be achieved using could eventually be achieved using continuous actuation. 

The power spectral density (Figure 8) shows that not much changes significantly as the continuous 
plasma control is applied. The growth of the turbulent kinetic energy (Figure 9) shows that there is 
upstream movement of the unsteady effects in the laminar shear layer that ultimately lead to the flow 
transitioning and reattaching as the geometry of the plasma actuation is increased. However, the primary 
growth rate of this turbulent kinetic energy remains constant for all of the cases, including the case of no 
plasma actuation. 





Figure 6. Time averaged velocity magnitude (z=0.1, pinching point, left) and Q-criterion (Q=100, 
right) for the cases of continual serpentine plasma actuation, varying the geometry of the actuator 
(from top to bottom — full, three quarters, half and one quarter). 
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Figure 7. Time and spanwise averaged coefficient of pressure, C,, on the surface of the airfoil (a) 
and time averaged local variations in the surface pressure coefficient - C,(x,z)-C,,.v, From top to 
bottom — Full, three quarters, half and one quarter (b-e). 


C. Pulsed actuation with the serpentine geometry — flow and surface features 

While the steady actuation shows limited degrees of success, when the actuation is pulsed at a 
frequency corresponding to instabilities in the shear layer, the reattachment is much more successful. The 
pulsed actuation causes the shear layer to transition and reattach. The boundary layer stays attached over 
the airfoil surface, the wake thins, and the drag experienced by the airfoil drops. 

The transitional vortical structures produced by the plasma actuation seem to be regular in structure 
These structures consist of a large hairpin structure and a more mobile smaller, hairpin structure 
component (Figures 10, 11). The large component is aligned with the spreading portion of the actuator 
geometry in the streamwise direction, while the smaller structure is aligned with the pinching portion of 
the actuator and lags behind the larger structure due to lower velocities closer to the airfoil surface. 


Power Spectral Density 
Power Spectral Density 





Figure 8. Power spectral density for the continuous actuation for x/c=0.5 (a) and 0.9 (b). 


These structures arise due to instabilities in the 
shear layer as it extends over the surface in the non- 
actuated portion of the duty cycle during the no plasma 
portion of the duty cycle. When they are initially 
formed, the larger structure (formed in the spreading 
NoControl region of the actuation) is very wide, and extends 


Full 


Three- Quarter upwards to the surface of the boundary layer, while the 


One-Quarter smaller structure (formed in the region of the pinching 
i oE Position °$ E actuation) is more compact and resides in the lower 
Figure 9. Growth of turbulent kinetic energy portion of the boundary layer. As these move 
along the surface of the airfoil with downstream, the smaller structure is thrust upwards 
continuous actuation. because of the small amount of streamwise vorticity 
Smak structure generated by the plasma actuator. This upward motion 
ng continues for the smaller vortex, which continues to 
grow as it convects downstream. As the two vortical 
structures pass over the trailing edge of the airfoil, the 
smaller structure has reached a similar height as the 
large vortical structure, if not reaching above it. 

It should be noted that while only the largest scale 
structures are easily visible for (Q=100), these 
structures are not single hairpin vortices, but packets of 
them, with additional hairpin vortices propagating off 
of the larger, initial structure as the packet moves 
downstream in a manner consistent with that of hairpin structures on a flat plate. For higher grid 
resolutions, the finer and finer structure associated with the vortices would become more visible and 
better resolved. 

From the comparison of the phase averaged and instantaneous flow, it can be seen that the two are 
very similar (Figure 13). The only noticeable difference is that the smallest unsteady structures have 
averaged out, most noticeable near the trailing edge and the wake. This suggests that the structures and 
their mixing are largely deterministic, and caused by the interaction of the shear layer and the plasma 
actuation. However, as these effects convect downstream, they transition into successively smaller and 
more random components. 

The phased averaged flow at the tail still exhibits very complex, yet seemingly periodic behavior 
(Figure 14). This indicated that a majority of the structures over the airfoil are deterministic in nature, and 
are repeated ever cycle. However, in the wake, it seems that there are more randomized components to the 
flow, and these components average out over a given number of samples. 
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Large structure 4 7 
Figure 10. Side view of the vortical structures 
from the pulsed actuation near the trailing 
edge. 


At the larger scale, these structures are periodic in 
the spanwise direction. This is due to the dictation of 
a dominant spanwise wave length by the serpentine 
plasma actuator. If this wavelength were varied, the 
length scale of the periodic behavior would change to 
match it. 

For the time averaged data, it can be seen that the 
flow does ultimately reattach, and in examining the 
Q-criterion, that some streamwise structures are 
created. These spanwise variations would not be 
present for a standard linear configuration (Figure 15, 
right). The large structures seen in the flow eventually 
average out, which exposes the underlying 
streamwise vorticity added by the plasma actuators 
(though some of this may be due to the regularity and 
length of the hairpin vortex legs). While the pattern of 
streamwise vorticies seen here seems very 
complicated, it should be recognized that the 
streamwise vortices farther in the leading portion of 
the chord leads to the hairpin vortical structures that Figure 11. Interaction of the hairpin vortices 
later dominate the flow. Additionally, since the from the phase averaged case with the full 
vortical structures aren’t directly connected to the electrode geometry during a complete cycle of 
actuator, the streamwise oriented structures are most Plasma actuation in increments of 7/2 starting 
likely a by-product of the actuation, a remnant of the from the top. The Q-criterion (Q=100) is 
disturbances introduced, not a direct contribution by Shown, colored by the velocity magnitude. 
the serpentine flow control. This does not mean that there is negligible addition of streamwise vorticity. 
The Q-criterion locates vortical structures based on the dominance of rotational stresses over irrotational 
stresses, not only the presence of streamwise vorticity, therefore, just as with the continual actuation, 
more shear than rotation is being added to the flow. 

The time averaged data also shows that there are slight regions of low momentum fluid over the 
surface of the airfoil (Figure 15, left). While these might be considered to be problematic, increasing the 
strength of the plasma actuator (either by additional force, or a denser, more compact force) would easily 
solve this. 

There also appear to be variations in the wake that vary with respect to the geometry of the plasma 
actuation. These variations are due to the size and shape of the hairpin vortices that are created. Two 
effects are occurring as the geometry of the device is changed. 

Examining the power spectral density (Figure 16) and of the resulting velocity fields within the 
boundary layer it can be seen, that behind the serpentine actuation at the mid-chord, the fluctuating kinetic 
energy available to the flow is moved out of the broadband spectrum and is largely put into the harmonics 
of the forcing frequency. Overall, the turbulent kinetic energy is large, but it is isolated into these 
harmonic frequencies. Closer to the trailing edge, the turbulent kinetic energy increases, but the 
magnitude of the turbulent kinetic energy remains smaller relative to the uncontrolled case. Furthermore, 
the energy stored in the peaks seems to decrease as the geometry of the actuator is straightened before 
increasing again for the flattest geometry tested here. 

Both of these effects show interesting behavior initially. At x/c=0.1, the strength of the disturbance 
increases as the serpentine geometry becomes more pronounced. However, all of the TKE and PSD(St=5) 
drop dramatically to roughly the same levels at x/c=0.2. Over the next half chord, the growth of the 
disturbance frequency is almost the same for all of the actuated cases, indicating that for the forcing 
distribution applied, the growth rate is roughly independent of the geometry of the plasma actuator. The 
growth of the TKE continues slightly farther down-stream, indicating that secondary effects are also 
coming into play, and that additional frequencies are also of importance. Examining the power spectral 
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densities shown previously, the first and second (St=10.0 and 15.0, respectively) are of greater or equal 
importance to the flow. 





Figure 12. Instantaneous velocity magnitude (z=0.1, pinching, left) and Q-criterion (right) for 
pulsed serpentine plasma actuation, varying the geometry of the actuator (from top to bottom - full, 
three quarters, half and quarter). 





Figure 13. Comparison of phase averaged (@=0) and instantaneous (bottom) velocity 
magnitude (left) and Q-criterion (Q=100, right) for the pulsed plasma actuation for the full 
geometry. 
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Figure 14. Phase averaged (@=0) velocity magnitude (z=0.1, left) and Q-criterion (Q=100, right) 
near the trailing edge of the airfoil 





Figure 15. Time averaged velocity magnitude (z=0.1, pinching point, left) and Q-criterion (Q=100, 
right) from above for the various geometries (from top to bottom - full, three quarters, half and 
quarter). 
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Figure 16. Power Spectral Densities for the pulsed actuation cases at (a) the midchord (x/c=0.5) and 
(b) near the trailing edge (x/c=0.9). 
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Figure 17. Growth of the turbulent kinetic energy (a) and disturbance frequency (b) along the 
surface of the airfoil. 
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Figure 18. Surface pressure (left), C, for the cases of pulsed serpentine actuation (a) and phase 
averaged local variations in the surface pressure coefficient - C,(x,zZ)-C) ay, at 6=0. From top to 
bottom — Full, three quarters, half and one quarter (b-e). 


D. Overall effects on lift and drag 

Though the variation of the geometry shows some impact on the lift and drag, it seems that the largest 
driver is whether or not the plasma forcing is pulsed or not, thus indicating that the unsteady stability is 
the most important thing to reattaching the flow (Table 3). Pulsing the actuation seems to bring a drop in 
Cı and Cp, with corresponding rise in <C,;/Cp> of roughly 19%, similar to the effects of reattachment for 
the linear actuator geometries'”''. Only marginal changes appear with continuously run plasma actuation, 
where <C;/Cp> increases only by 1.0%. 

In examining the transient forces experienced by the airfoil under steady actuation, it can be seen that 
there is at least one important frequency present in the forces (Figure 19). This frequency is slightly lower 
than the St=.05 frequency seen with no actuation, suggesting that the shedding frequency of (now three 
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dimensional) boundary layer has changed. Furthermore, a slower level of unsteadiness seems to still exist 
for the continuously actuated cases. Again, this may suggest that some slow behavior exists that does not 
exist for the un-actuated or the pulsed actuation simulations. 

The lift and drag for the cases of pulsed actuation also shows oscillatory behavior, but that behavior is 
consistent with the St=5 forcing frequency (Figure 20). Additionally, aside from a slight increase for the 
fullest serpentine geometry, the magnitude of the oscillations increases as the geometry compresses, 
possibly resulting from the more coherent vortical structures in the boundary layer passing over the 
trailing edge all at once (similar effects with spanwise vortical structures were seen previous work'””’). 
These oscillations are also smaller than those of the uncontrolled or of the continuously actuated case. 


Case <C > <Cp> <C> <C,/Cp> Cr’ Cp í Cyr í 

Baseline 0.58926 0.0246 -0.0485 23.998 0.0087 0.00091 0.0042 
Full — Pulsed 0.54782 0.0193 -0.0303 28.433 0.0029 0.00028 0.0015 
Full - Cont. 0.56207 0.0234 -0.0378 23.986 0.0079 0.00076 0.0032 
Three Quarters — Cont. 0.56061 0.0229 -0.0372 24.508 0.0058 0.00058 0.0026 
Half — Cont. 0.57124 0.0235 -0.0410 24.297 0.0064 0.00065 0.0030 
One Quarter — Cont. 0.56827 0.0235 -0.0409 24.190 0.0110 0.00094 0.0043 
Three Quarters — Pulsed 0.54551 0.0191 -0.0304 28.530 0.0018 0.00022 0.0010 
Half — Pulsed 0.55076 0.0192 -0.0310 28.476 0.0030 0.00035 0.0020 
One Quarter — Pulsed 0.55358 0.0193 -0.0318 28.678 0.0036 0.00044 0.0023 


Table 3. Average (<...>) and standard deviations (‘) of the airfoil forces for the various cases run 
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Figure 20. Time resolved C, (a) and Cp (b) with pulsed plasma actuation. 
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V. Conclusions 


Simulations were performed in order to better understand the effects of geometric variation using 
plasma actuation. While these simulations achieved the desired result of reattaching the flow and show 
that the fullness of the serpentine geometry impacts the flow, but that other effects are likely more 
important. 

While serpentine geometry plasma actuation shows promise as a tool for flow control under other 
circumstances, it appears that for the forcing distibution tested, the geometric curvature of the electrodes 
does not provide additional control authority of separation mitigation. For this particular case, the stability 
of the separated shear layer seems to be the most important parameter, and pulsing the flow ultimately has 
the largest effect on reattachment. 

While this device does not show significantly more or less control authority as compared to the 
traditional linear plasma actuator, it did show a slightly different method of control authority because it is 
able to apply a characteristic wave number in the spanwise direction, and influence the transition and 
reattachment in that manner. Future efforts could be directed towards replacing spanwise distributed 
roughness or some other flow control method. Further exploration of the effects of expanding into these 
novel geometries under canonical flow conditions may help to identify conditions where these devices 
may be the optimal flow control device. 
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